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Recently evidence has emerged for enormous features in the 7-ray sky observed by the Fermi- 
LAT instrument: bilateral 'bubbles' of emission centered on the core of the Galaxy and extending 
to around ±10 kpc above and below the Galactic plane. These structures are coincident with a 
non-thermal microwave 'haze' found in WMAP data and an extended region of X-ray emission 
detected by ROSAT. The bubbles' 7-ray emission is characterised by a hard and relatively uniform 
spectrum, relatively uniform intensity, and an overall luminosity ~4x 10 37 erg/s, around one order 
of magnitude larger than their microwave luminosity while more than order of magnitude less than 
their X-ray luminosity. Here we show that the bubbles are naturally explained as due to a population 
of relic cosmic ray protons and heavier ions injected by processes associated with extremely long 
timescale ( > 8 Gyr) and high areal density star-formation in the Galactic center. 
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A recent analysis [TJH] of Fermi-LAT [3] 7-ray data has 
revealed two, enormous bubble-like emission features cen- 
tered on the core of the Galaxy and extending to around 
±10 kpc above and below the Galactic plane. At lower 
Galactic latitudes these structures are coincident with 
a non-thermal microwave 'haze' found in WMAP 20-60 
GHz data [H [5] and an extended region of diffuse X-ray 
emission detected by ROSAT [BJ. 

A natural explanation of these structures [2] would 
be that they are due to the same population of highly- 
relativistic ( > 50 GeV) cosmic ray (CR) electrons which 
synchrotron radiate at multi-GHz frequencies and simul- 
taneously produce > 1 GeV 7-rays through the inverse 
Compton (IC) process. However, given the severe ra- 
diative energy losses experienced by electrons, the hard 
spectrum, uniform intensity, vast extension, and energet- 
ics of the bubbles render the origin of this particle popu- 
lation extremely mysterious (see fig. [T]) [fl El SI El E] . In 
particular, transport of >TeV, IC-radiating electrons to 
the requisite distances from the plane would require ve- 
locities of > 0.03 c, too fast for a Galactic wind (this is a 
conservative lower limit as only electron energy losses on 
the 2.7 K CMB arc accounted for). If diffusive, a diffusion 
coefficient of ~ 10 31 cm 2 /s would be required for E e = 1 
TeV, 1-2 orders of magnitude larger than the Galactic 
plane value. One might postulate an in-situ electron ac- 
celeration/injection process to surmount these difficulties 
PQ but here spectral considerations present a severe test. 
The ~ E~ 7-ray spectrum might be due to IC emission 
from a cooled ~ E~ 3 electron population but there is 
a robustly-detected |l a hardening in the 7-ray spectrum 
below ~ 1 GeV. In order to produce such a break either 
a unique (over the age of the Galaxy) injection event of 
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age ~ 10 6 years or a sharp, ~1 TeV low-energy hard- 
ening or cut-off in the injection spectrum of electrons is 
required. The former seems unlikely as there are no in- 
dications of such an event occurring over this timescale 
in the Galactic center (GC). The latter cannot be ex- 
cluded in principle: an in-situ acceleration process (in- 
volving, e.g., stochastic acceleration or magnetic recon- 
nection) that produced a spectrum harder than ~ E^ 1 
for E e <1 TeV would suffice, though such seems ad hoc. 

In contrast to the difficulties presented by the elec- 
tron/IC case, we show below that a CR proton popu- 
lation - associated with extremely long timescale star- 
formation (SF) in the GC and injected into the bub- 
bles by a wind- can naturally explain the 7-ray struc- 
tures provided the protons are trapped for timescales ap- 
proaching 10 10 years. Our basic picture is that, granted 
this timescale, the majority of the power going into non- 
thermal protons is lost into pp collisions on the bubbles' 
low-density plasma (itself the 'wind-fluid' also injected by 
GC SF) and subsequently reprocessed into 7-rays, elec- 
trons, positrons and neutrinos. Because the protons are 
trapped in our scenario their steady-state distribution 
mirrors their injection spectrum. Requiring that the bub- 
bles' hard-spectrum and overall (1-100 GeV) luminosity 
be reproduced implies a hard spectrum population of CR 
protons of (time-averaged) power ~ 10 39 erg/s injected 
over multi-Gyr timescales. Given that our independent 
[TO] work suggests that i) the GC launches > 10 39 erg/s 
in hard-spectrum CRs on a strong outflow, ii) the mor- 
phology of the bubbles clearly privileges the GC [T] , and 
iii) the GC is perhaps the single, spatially-localized site 
in the Galaxy where SF over multi-Gyr timescales is as- 
sured 8, 9 J this SF ultimately offers a compelling expla- 
nation of the 7-ray and - incidently - microwave and 
X-ray phenomenology of the bubbles as we now explore. 

Cosmic ray hadrons (in principle protons and heavier 
ions, but, for simplicity, we refer to the dominant protons 
below) undergo collisions with ambient matter creating 
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FIG. 1: Timescales for CR protons and electrons in 
the bubble medium. Key:- 't p riax, ■. maximum cooling 
time, l/(d In E/di)), for protons given ionization and pp col- 
lisions in the ur — 0.005 cm -3 bubble plasma. The vertical 
solid lines show the approximate mean proton primary energy 
for daughter 7-rays of 1 GeV and 100 GeV. 't™?™': minimum 
formation time (~ 8 Gyr) for the bubbles in order that the 
system reach saturation. ^t™ ax ': maximum timescale (adopt- 
ing Ub <C Ucmb or B < 3 x 10~ 6 G) of electron cooling 
due to the combination of ionization and bremsstrahlung in 
the plasma and IC-scattering on the CMB. The solid hori- 
zontal line shows the timescale for transport of electrons out 
to the full extent of the bubbles on a (very fast) 1000 km/s 
wind. The vertical dashed lines show the energies required for 
an electron to IC scatter a CMB photon to 1 and 100 GeV. 
Much higher energy electrons than protons are required to 
produce a daughter 7-ray of the same energy. 



daughter mesons (mostly pions), the neutral component 
of which decays into 7-rays. This process explains most 
of the diffuse 7-ray emission detected at > 100 MeV en- 
ergies from the Galactic plane. The hadronic 7-ray lu- 
minosity of a region scales linearly in its thermal (gas) 
and relativistic (CR hadron) populations, thus, assum- 
ing that high energy CRs are highly penetrating, it is 
generally taken that astrophysical markers of ambient 
gas density should trace hadronic 7-rays. This heuris- 
tic does not always apply, however. If the timescale 
for particle acceleration events in a system is smaller 
than all other (energy loss or escape) timescales the sys- 
tem will be in (quasi) steady-state. If the timescale for 
particle energy loss (via hadronic collisions here) is less 
than the escape time the system qualifies as a 'thick tar- 
get'. If both conditions pertain the system is in satu- 
ration. The (hadronic) 7-ray luminosity of a region is 
~ Np/tpp^K ! where the N p is the region's steady- 



state proton population and tpp^^o is the timescale for 
neutral pion production in such collisions. In satura- 
tion, however, N p ~ Q p t pp (because energy loss through 
pp collisions is the dominant loss process) and - reflect- 
ing the almost equal production of 7r°,7r + , and 7r~ in 
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FIG. 2: Fermi bubble 7-ray spectrum [T]. The error 
bars are la. The solid curve shows the best-fit spectrum 
(oc E~ 2A ), the dotted curve shows the steepest reasonable 
spectrum (oc Ep 23 ), and, for comparison, the dashed curve 
is the spectrum expected were the bubbles suffused with pro- 
tons having the oc E~ 2 ' 7 distribution of the Galactic disk. 



L«i ~ Qp/i- Thus, in saturation about a third of the 
power injected into relativistic CRs will emerge in j-rays 
(of all energies), independent of in-situ gas density, in- 
teraction volume, and CR injection time. 

Assuming saturation, a hadronic scenario reproduces a 
number of the aspects of the phenomenology of the Fermi 
bubbles: i) the reported [UH] hard 7-ray spectrum is ex- 
plained; consider the the situation in the Galactic plane. 
Here diffusive confinement of the CRs leads to a steep- 
ening of the steady state spectrum to oc E~ 2 - 7 . In stark 
contrast, there is no energy- dependent confinement effect 
in the bubbles. So, given the almost energy-independent 
pp loss time, we see the spectrum of the CRs as injected at 
their acceleration sites (evidently oc E~ 2A ) mirrored by 
the bubble 7-rays. ii) Decay kinematics, however, do im- 
ply a bump at low energy in the 7-ray spectrum around 
E~ t = m n o/2 ~ 70 MeV. This is mapped to a down- 
turn below ~GeV on a spectral energy distribution plot 
(see fig-[2j n.b. the lowest and highest energy data points 
suffer from large systematic uncertainties [1] in the sub- 
traction of back/foreground emission); this explains the 
afore-mentioned down-break in the bubble 7-ray data, 
iii) In saturation circumstances conspire to establish a 
constant volume emissivity - regions of higher gas den- 
sity necessarily have attenuated steady-state proton pop- 
ulations because of extra cooling by the gas. This effect 
tends to give a uniform intensity of emission, as observed, 
despite the expected variations in target gas density. 

To satisfy the requirement that the system be in satu- 
ration we require that the timescale associated with en- 
ergy loss through pp collisions (the dominant loss process 
at relevant energies) is less than the timescale over which 
CRs have been injected: 
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Given the inelasticity satisfies k pp ~ 0.5, a pp ~ 5 x 10~ 26 
cm 2 , and an plasma density ~ 0.01 cm -3 1], this condi- 
tion implies t in j > 5 x 10 9 (n#/0.01 cm -3 )" 1 yr. 

It is remarkable that the rate of SF in today's GC can 
supply the power and mass required to explain the bub- 
bles' 7-ray emission. IRAS data [12] imply that the inner 
1.5° (in diameter) around the GC emits a total infrared 
luminosity Ltir — 1.6 x 10 42 erg/s which implies [T3] a 
star formation rate for the region of SFR ~ 0.08 M /yr. 
Such a SFR 7 in turn, implies [14) a supernova (SN) rate 
throughout the same region of 0.04/century (uncertain 
by ~2) leading to a total power injected by supernovae 
of E SN ~ 1.3 x 10 40 E 51 erg/s, where E 51 = E SN /1Q 51 
erg is the mechanical energy released per SN. Assuming 
[T5"[ ~ 10% of a SN's mechanical energy goes into accel- 
erating non-thermal particles, this implies GC generates 
a power of ~ 1.3 x 10 39 erg/s in CRs. 

From Fermi data [16] we infer that the 1° x 1° field 
around the GC emits ~ 3 x 10 36 erg/s s in > 1 GeV 
7-rays, much less than the > 10 38 erg/s expected were 
the system calorimetric and given its SN rate (the GC 
GeV emission is substantially polluted by line-of-sight 
and point source emission in any case; the 7-ray flux at 
~ TeV energies reported by HESS [T7] also indicates the 
system is far from calorimetric [ID]). The GC, then, loses 
~ 10 39 erg/s in hard-spectrum CRs into the Galactic en- 
vironment. In comparison, the bubbles emit 4 x 10 
erg/s in 1 to 100 GeV 7-rays [T], implying that CR pro- 
tons in the energy range 10 to 1000 GeV lose ~ 1.2 x 10 38 
erg/s in pp collisions or ~ 3.6 x 10 38 erg/s integrating the 
~ E~ 2 distribution from GeV to 10 6 GeV. Accounting 
for ionization losses by sub-relativistic protons and adia- 
batic energy losses at all energies, bubble protons lose a 
total ~ 10 39 erg/s in steady state, precisely accounting 
for the CR power injected at the GC. 

A prime candidate for a process that removes most 
of the GC's CRs is a 'super-wind' driven by the same 
SF processes ultimately responsible for the CR acceler- 
ation. Such winds are detected emerging from the nu- 
clei of many star-forming galaxies [IB] . A detailed ac- 
counting |11] - taking into account mass loss from stars 
and injected by supernovae - finds that the mass in- 
jected into this wind by the inner 1.5° region is M wind — 
(0.02 — 0.03) M Q /year. The asymptotic speed of such a 

wind will scale as v winc i ~ J 2 Esn /M win( j which evalu- 
ates to 1200 km/s (for an adiabatic outflow), larger than 
the gravitational escape speed from the region (~ 1000 
km/s [19 ). There are, however, both strong empirical 
[2"T] and theoretical [5D] indications that, given radiative 
losses, such an outflow should stall at a height < 15 kpc, 
consistent with the ~10 kpc half-height of the bubbles. 

With the above mass injection rate, we can re- visit the 
the minimum formation timescale for the bubbles under 
the assumption that the CR target gas is precisely that 
carried out of the GC on the putative wind. The gas 
density in the bubbles is Hh — Mbbi x Unj/Vbhl where 
-^bbi™ < Mai < -^wind, given that the net mass growth 



of the bubbles has to allow for plasma cooling and falling 
back to the plane (M^f 1 = 0.004 M /yr corresponds to 
tinj = 13 Gyr). In concert with eq. [I]we find 
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This gas is injected at the base of the wind dominantly 
in plasma form and will radiate an amount of power in 
thermal bremsstrahlung X-rays controlled by the plasma 
density (already determined) and temperature of 
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where A(T) is the plasma cooling coefficient [55] and the 
last inequality follows from the preceding equations. In 
steady state, conservation of energy demands that the 
plasma luminosity be less than the power injected at the 
base of the wind. This consideration - together with 
ROSAT observations of the region JH] [23] (which reveal 
a diffuse, < 10 7 K plasma coincident with the bubbles 
[TJ H] at relatively lower Galactic latitudes) - imply that 
the plasma is extremely hot (as previously contemplated 
PQ): > 3x 10 6 K and of density (0.004-0.006) cm-' 3 . This 
temperature scale is similar to that observed for 'super- 
winds' from the nuclei of star- forming galaxies [23] . 

Two effects may operate to ensure the plasma's high 
temperature. Firstly, it is injected at very high tempera- 
ture (10 K), implying long radiative cooling times, ~ 
Gyr (there have been persistent, though disputed [2"5] . 
claims that observations with Chandra [15] , SUZAKU 
[26] . and previous X-ray instruments [27 suggest a 
(6 — 9) x 10 7 K, diffuse plasma in the ~ 100 pc around 
the GC). Secondly, the plasma may be (re-)heated by 
thermalization of initial bulk motion: the kinetic energy 
of each thermal proton in an 800 km/s outflow corre- 
sponds to a temperature of 2 x 10 7 K. Moreover, given 
the long formation timescales, coloumbic processes (act- 
ing on timescales <C 10 7 yr) may distribute this final 
thermal energy between plasma protons and electrons. 

A hadronic 7-ray scenario predicts the injection of rel- 
ativisitic secondary electrons (and positrons) throughout 
the bubbles, synchrotron radiating on ambient magnetic 
fields. In fact, in the calorimetric limit, the kinematics 
of the charged and neutral pion decay chains (and the 
oc E 2 dependence of the synchrotron critical frequency) , 
imply that vL v {synch) is a significant fraction - up to a 
quarter (depending on synchrotron vs. non-synchrotron 
energy loss rates) - of vL v { / k°). For the best-fit 7-ray 
spectrum (fig. [2| we find a 20-60 GHz secondary electron 
synchrotron luminosity of <~ 2 x 10 -36 erg/s for B > 10~ 5 
G (implying near energy density equipartition between 
the bubbles' magnetic field, CRs, and plasma). This is 
consistent with observations: the luminosity of the haze 
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over the same WMAP bands (20-60 GHz) is (1-5) x 10 36 
erg/s [4] to be contrasted with the ~ 4 x 10 37 erg/s in 7- 
rays [T] . The energy-independent transport timescale as- 
sociated with the putative wind (cf. [53]) can explain the 
hard spectrum of the microwave haze (where it is mea- 
sured inward of ~ 18° [7]). If it transpires that the mi- 
crowave spectrum is hard out to the full, ~ 10 kpc extent 
of the 7-ray bubbles this would be unexplained within 
our scenario and imply an additional electron population 
(injected or accelerated in situ as previously discussed). 

Our scenario requires CR proton trapping for multi- 
Gyr timeframes. Diffusive confinement could achieve this 
alone for a diffusion coefficient 1-2 orders of magnitude 
smaller than the Galactic plane value (cf. the electron/IC 
scenario). This is reasonable: X-ray observations [50] 
show that the plasma of the Galactic bulge is extremely 
turbulent. Moreover, a search for polarized emission in 
the WMAP data [31] has yielded a negative result, rec- 
oncilable with this emission's being due to synchrotron if 
the magnetic field structure is highly tangled [7]. If the 
claimed [1] sharpness of the bubble edges is correct, it 
is likely that a magnetic draping effect [25J around the 
super-Alfvenic GC wind also helps trap the protons as 
apparently occurs in fossil cluster radio bubbles [3"2"] . 

The GC wind probably also has an important role in 
transporting positrons from the Galactic core out into the 



Galactic bulge which is traced by the electron-positron 
annihilation line emission |33j . This naturally predicts 
511 keV emission out to ~ 1 kpc from the GC, as ob- 
served, despite the fact that positrons must be injected 
into the ISM at energies of only a few MeV [3"4ll3"5] . 

Finally, we predict there should be extended, ~ TeV 7- 
radiation surrounding the Galactic nucleus on similar size 
scales to the bubbles with an intensity < E^-F 7 (TeV) ~ 
10~ 9 TeV cm -1 s _1 sr _1 which should make an interest- 
ing target for future 7-ray studies. Likewise, the region 
is a promising source for a future, Northern Hemisphere, 
km 3 neutrino telescope for which we estimate (assuming 
a 7 = 2.0 proton spectrum, cut-off at 1 PeV) ~40 signal 
events above 10 TeV per annum (vs. ^100 background 
events implying a 4<r detection in one year). Of course, 
steepening or a lower-energy cut-off in the intrinsic pro- 
ton spectrum or, indeed, loss of confinement at higher 
proton energies would reduce the significance of the bub- 
bles as either TeV 7-ray or neutrino sources. 
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